Ising Spin Glasses : interaction distribution dependence of the critical exponents 
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Continuous phase transitions are catalogued into universality classes, families of systems having 
identical values of all the exponents governing the critical behaviour of their different physical 
properties. Numerical simulations have been carried out on Ising Spin Glasses in dimension three, 
using a technique where corrections to finite size scaling can be controled. The data show that 
the critical exponents vary strongly as a function of the kurtosis of the interaction distribution, a 
parameter which from the standard point of view should not be pertinent. This observation implies 
that for spin glasses the renormalization group analysis should not be approached in a the same way 
as in the case of canonical second order transitions; a much richer structure of universality classes 
would appear to exist for spin glasses. 
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The renormalisation group (RG) provides an explana- 
tion of the physical origin of the critical exponents and of 
the universality classes for standard second order transi- 
tions which is one of the most remarkable achievements 
of statistical physics. The university rules state that the 
critical exponents depend only on a small number of basic 
parameters, typically the dimension of space d and the 
number of order parameter components n |jj . Close to a 
transition, large fluctuating correlated clusters form; uni- 
versality can be viewed as the consequence of the large 
scale fractal structure of these clusters being indepen- 
dent of the details at the microscopic level of the inter- 
actions between the elements (spins, atoms,...) making 
up the system. The few known exceptions to universality 
concern mostly rare but well understood marginal cases, 
such as certain regularly frustrated spin systems in two 
dimensions, where critical exponents vary continuously 
with the value of a control parameter. 

Glass transitions remain enigmatic despite a consid- 
erable theoretical and experimental effort S; even the 
existence or otherwise of a thermodynamic transition in 
structural glasses is still the subject of intense debate. 
Spin Glasses can be considered the magnetic analogues 
of the familiar structural glasses; the canonical examples 
are made up of local magnetic moments which interact 
randomly with each other. For technical and conceptual 
reasons the spin glass transition is much easier to study 
than structural glass freezing; in particular powerful nu- 
merical techniques can be readily applied for Ising Spin 
Glasses (ISGs). Although spin glass freezing is an un- 
orthodox form of transition as compared with standard 
second order transitions, extensive experimental and nu- 
merical work has demonstrated that these freezing tran- 
sitions show true continuous transition critical behaviour 
albeit with exponents which are very different from stan- 



dard second order transition values j|, ^ §|. Deter- 
mining transition temperatures and critical exponents to 
high precision at spin glass transitions, either experimen- 
tally or through simulations, is however a notoriously dif- 
ficult task. We estimate values for the freezing temper- 
atures and critical exponents of ISGs in dimension three 
using numerical techniques where systematic errors can 
be kept carefully under control. Confirming earlier work 
0, H , we find sets of critical exponent values which vary 
strongly with the form of the interaction distribution, in 
contradiction to the conventional universality rules. 

Problems that must be faced for numerical simulations 
in ISGs include the agonisingly slow equilibration rates, 
and the need to average over numerous microscopically 
inequivalent samples to obtain a result representative of 
the true mean over all possible samples with a given ran- 
dom interaction rule. Calculations have necessarily been 
restricted to moderate size samples, finite size scaling 
rules being used to evaluate the thermodynamic limit 
ordering temperature T g and the associated critical ex- 
ponents |^] . An important caveat concerns deviations 
from finite size scaling which can corrupt the values of 
the critical parameters estimated from data taken only 
over a narrow range of small to moderate sizes [To) . 

In Ising Spin Glasses (ISGs) a basic property which 
is measured numerically is the mean autocorrelation or 
memory function averaged over spins Si as a function of 
time t : 



q(t) =< Si(t)Si(0) > 



(1) 



The freezing temperature T g is defined as the tempera- 
ture below which q(t) remains non-zero to infinite t in the 
thermodynamic limit Jl~T[ ] . It can be estimated from the 
divergence of the relaxation time for very large samples 
Q. Alternatively, the standard finite size scaling tech- 
nique for determining T g from measurments on moderate 
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size samples is through the dimensionless Binder param- 
eter ratio of the moments of the equilibrium distribution 
ofg(f) @: 

g L = 1/2 [3- < q A > I < q 2 > 2 ] (2) 

The curves (T) for different sizes L should all intersect 
at T g . Unfortunately in 3d ISGs intersections need very 
high quality statistics to be visible at all ^2) and can 
be subject to severe deviations from finite size scaling 

H & 

An alternative approach which has been used to iden- 
tify T g and associated critical exponents Q| involves the 
combination of measurements of three parameters. 

First, equilibrium measurements of the spin glass sus- 
ceptibility 

XSG = L d [<q 2 >] (3) 

on samples of different sizes L obey the scaling rule 

Xsg/L 2 ~ L- 71 (l + const *L- U + ) (4) 

at T g with n being the standard static critical exponent 
Q and uj being the leading static correction to finite size 
scaling exponent. Below T g it turns out that the same 
functional form can be taken to hold with a temperature 
dependent effective value r](T). In contrast to the case of 
the two parameters that we will discuss next, direct xsg 
measurements require long preparatory anneals or the 
use of parallel tempering methods to assure true equilib- 
rium and so are much more demanding numerically than 
measurements on non-equilibrium parameters. The need 
to attain true thermal equilibrium in finite computing 
time means that the maximum sample size for the sim- 
ulations is in practice limited to L ~ 20. Fits to high 
quality data (including results at small L) and allowing 
for the correction to scaling factor can nevertheless pro- 
vide accurate values for r](T) together with an estimate 
for the exponent w ||. 

Secondly, following a preliminary anneal at a temper- 
ature T near T g over a long waiting time t w , the initial 
relaxation of the autocorrelation function takes the form 

q(t-t w )=X(t-t w )- x ^ (5) 

As long as the anneal time t w is much longer than the 
measuring run time t — t w , the value of x(T) observed is 
equal to the equilibrium value (corresponding to infinite 
t w ), and the value x(T) which can readily be obtained on 
reasonably large sized samples is essentially equal to the 
infinite size limit [TJJ] . Direct checks have been made of 
the independence of the measured x(T) on sample size, 
and the absence of corrections to finite size scaling be- 
yond a few Monte Carlo Steps per spin (MCS) Q . Hence 
precise values of x{T), for all intents and purposes equal 
to the limiting equilibrium values in the thermodynamic 



limit, can be obtained by averaging over data on a suffi- 
cient number of large samples without the need to anneal 
to true equilibrium. At T = T g , x(T g ) is related to the 
standard static and dynamic critical exponents rj and z 
through |sj 

x{T g ) = {d-2 + r 1 )/2z (6) 

Finally, the out of equilibrium time dependent spin 
glass susceptibility for a sample quenched from infinite 
temperature and annealed towards equilibrium at tem- 
perature T is defined by 

Xsg(T) = [< S?(t)S b i(t) > 2 ] (7) 

where a and b are two independent replicas of the same 
sample initially at infinite temperature at t = 0. Then 
for times greater than some microscopic time 

X ' SG (t) = At h (T)(l + B{r w/z + higherorder) (8) 

where h(T) is the "Huse exponent" |l5j and w is the lead- 
ing dynamic correction to scaling exponent [ pl| (which in 
principle does not have to be equal to the leading static 
correction exponent w Q). The corrections at short times 
are the analogue of finite size corrections, as the clusters 
of correlated spins that are gradually building up with 
increasing annealing time have a t dependent finite size, 
even if the sample size L can treated as effectively infi- 
nite. The factors A and B may vary with T. At T = T g 
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h(T g ) = (2 - T,)/z (9) 

Here, by definition there is no equilibration anneal before 
a measuring run to determine h(T), and large samples 
can readily be used, so it is relatively easy to ensure that 
for the time range used the system is far from the satura- 
tion limit for that particular size, see jL7|. (Direct tests 
can be made by comparing data taken for different large 
values of L; when the h(t) values become L independent 
then one is in the infinite size limit). By averaging over a 
sufficient number of runs on independent samples, accu- 
rate results for of < Xsg(^) > can ^ e obtained, so h{T) 
can be estimated to good precision from a fit including 
the leading correction to scaling. Figure 1 shows exam- 
ples of a plot of data with a fit for the binomial interac- 
tion distribution. It turns out that the correction can be 
measured but is much weaker in the other two cases. 

Suppose we now make measurements of the three in- 
dependent parameters x(T{), h(Ti) and rj(Ti) at a series 
of test temperatures Ti around T g . We can define a new 
effective parameter h*(Ti) through 

h*{Ti) = 2x{Ti){2 - fj(Ti))/(l + rj(T,)) (10) 

and we can plot both h(Ti) measured directly and the 
derived h*(Ti) against x(Tj), with Tj as an implicit pa- 
rameter. 
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FIG. 1: The non-equilibrium spin glass susceptibility Xsg(^) 
as a function of time t for the binomial ISG at temperature 
T — 1.2. The fit is to equation 8 with the Huse exponent 
h(T) — 0.398 and the leading correction exponent w/z = 0.61. 



For one given system, because of the definitions of h 
and h* at T g in equations 6, 9 and 10, consistency dictates 
that the curves for h(x) and for h*(x) must intersect at a 
point which corresponds to Tj = T g and which therefore 
lies at the point x — x(T g ) and h = h* = h{T g ). If stan- 
dard universality holds, the values of the static critical 
exponents Tj(T g ) and of the dynamic critical exponents 
z{T g ) at the appropriate critical temperatures should be 
identical for all short range interaction ISGs in a given 
dimension. The form of the random interaction distribu- 
tion should not be pertinent. Therefore the curves h{x) 
and h*{x) for all ISGs in a particular dimension should 
intersect at one single point corresponding to the univer- 
S <d[h(T g ),x(T g )}. 

We have made measurements on standard three dimen- 
sional ISGs on simple cubic lattices with toroidal bound- 
ary conditions. The three systems chosen have random 
near neighbour interaction distributions which are bino- 
mial, Gaussian, and decreasing exponential respectively. 
The data for x(T) and h(T) are from our own simula- 
tions. To obtain these parameters samples of maximum 
size L — 28 were studied, with an average being taken 
over 6600 independent samples. For the x{T) measure- 
ments preliminary anneals were performed to at least 10 6 
MCS, and the run times used were up to 10 4 MCS. For 
the binomial interaction distribution and for the Gaus- 
sian interaction distribution, x(T) values are in excellent 
agreement with data shown in Q and [[u] at all tem- 
peratures where direct comparisons can be made. There 
do not appear to be published measurements for the de- 
creasing exponential distribution. For the h(T) measure- 
ments, runs were made to 10 4 MCS. Our binomial in- 
teraction data are in good agreement with but are more 
accurate than those of |l5|, [l8|]. We are not aware of 
published h(T) data on the other systems. For the rj(T) 



values we rely on accurate size dependent spin glass sus- 
ceptibilty \sg results by ||, [ll| and ||(J for the binomial 
distribution up to L — 24 and the Gaussian distribution 
up to L = 16 respectively. The ~q(T) data for the de- 
creasing exponential distribution correspond to our own 
measurements of x(^) U P to only L = 8. The statistical 
precision of the x(T) and h(T) values is about ±0.005 for 
all three systems, and as we have explained the residual 
systematic error due to incorrect extrapolation to infi- 
nite size can be taken to be small. The precision of the 
rj(T) values is about ±0.02 for the binomial and Gaussian 
interaction systems, with the uncertainty being mainly 
statistical with a possible contribution from incomplete 
equilibration. The uncertainty is greater for the decreas- 
ing exponential case because of relaxation rates for equi- 
libration which are even slower than for the other inter- 
actions, limiting the sizes L which we could anneal to full 
equilibrium. 

The results are presented in Figure 2. For each of the 
three systems individually, the critical values x(T g ) and 
h(T g ) can be estimated directly from the [h(x),h*(x)] 
intersection point. The corresponding critical values r\ 
and z follow immediately from equations 4 and 6. As 
we know x(T), the values of the critical temperatures T g 
follow also from the positions of the intersection points. 
The values estimated are given in Table 1. The values 
of the critical temperatures can be compared to other 
published estimates §, |, g, [l0|, @ |o). 
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FIG. 2: The effective exponents h(x) (closed symbols) and 
h*(x) (open symbols) as functions of the exponent x, with 
the temperature T as an implicit parameter. Circles : bino- 
mial distribution, triangles : Gaussian distribution, squares 
: decreasing exponential distribution. The numbers indicate 
temperatures for the different distributions. The intersection 
points of h(x) and h*(x) correspond to the critical tempera- 
ture T g for each distribution. 

A cursory glance at Figure 2 shows that the intersec- 
tions of the h(x) and h*(x) curves for the three different 
systems do not coincide at one single point in the [h, x] 
plane, so the naive universality rule, which would have 
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the ISG transitions lying within a single universality class 
whatever the form of the interaction distribution, is vio- 
lated. 

If we consider the tabulated values, it can be seen that 
the various exponents do not vary in a random way. Go- 
ing from one system to the next, as the kurtosis of the 
interaction distribution increases and the T g values di- 
minish, r] becomes more negative and z increases. This 
is just the tendency one would expect if the exponents 
were tending regularly towards the T g = values, which 
are z = oo, and r\ = — 1 in dimension 3. We will discuss 
elsewhere the behaviour of the correction to scaling ex- 
ponent w |2ll , but we can note that the analyses indicate 
high values of w. In [0 a high cu estimate was linked to 
the fact that dimension 3 is far below the upper critical 
dimension, 6. High values of w and oj are favourable, as 
they mean that corrections to finite size scaling die out 
relatively rapidly with t and L respectively. 



TABLE I: Estimates of T g , 77 , 2 and w/z for the 3d ISG with 
different interaction distributions 
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Interaction 


T s 
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z 


w/z 


Binomial 


1.190(15) 


-0.20(2) 


5.85(10) 


0.48(3) 


Gaussian 


0.920(15) 


-0.42(2) 


6.45(10) 


0.62(3) 


Decreasing Exponential 


0.70(5) 


-0.54(5) 


8.8(3) 


0.95(2) 



Universality has a strong aesthetic appeal, and the 
general physical argument outlined in the introduction 
would seem very robust ; nevertheless in a different con- 
text there are examples of systems with many attractors 
where critical exponents vary continuously as a function 
of a control parameter [^2| . Evidence has been found for 
universality in Migdal-Kadanoff (MK) spin glasses |p3| ; 
however MK spin glasses are known to have a "droplet" 
structure p4| , while there are strong arguments in favour 
of a more complex structure in finite dimension ISGs (see 
for instance [^5| ). In the finite dimension spin glass con- 
text the e expansion below the upper critical dimension 
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in ISGs has proved totally intractable 26 
p7| state that the classical tools of RG analysis are not 
suitable for spin glasses. They introduce a non-standard 
coarse graining technique that is performed on the over- 
lap probability measure. Within this generic approach 
a much richer structure of universality classes could ap- 
jrj spin glasses, and in particular the form of the 
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interaction distribution could play a fundamental role. 



Thus in ISGs there could well exist fundamental reasons 
to expect changes of critical exponents with interaction 
distributions like those we observe empirically. 
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